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Abstract 

We point out the possibility that reactor measurement of ^13, when combined with high-statistics 
appearance accelerator experiments, can detect leptonic CP violation. Our proposal is based on 
a careful statistical analysis under reasonable assumptions on systematic errors, assuming 2 years 
running of the neutrino mode J-PARC Hyper-Kamiokande experiment and a few years running 
of a reactor experiment with lOOton detectors at the Kashiwazaki-Kariwa nuclear power plant. We 
show that the method can be arranged to be insensitive to the intrinsic parameter degeneracy but 
is affected by the one due to unknown sign of Am-g^. 
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I. INTRODUCTION 



After the pioneering and the long-term extensive efforts in the atmospheric the solar 
0], the accelerator and the reactor j^] experiments, we have grasped the structure of 
lepton flavor mixing in the (2-3) and the (1-2) sectors of the Maki-Nakagawa-Sakata (MNS) 
matrix j^. Now we are left with the unique unknown (1-3) sector of the MNS matrix, 
in which there live the third mixing angle ^13, which is known to be small 0], and the 
completely unknown CP- violating leptonic Kobayashi-Maskawa phase S 01 • 

Detecting leptonic CP violation is one of the most challenging goals in particle physics. A 
popular method for measuring the CP-violating phase is by long-baseline (LBL) accelerator 
neutrino experiments using either conventional neutrino superbeam [sl. [lolllll| . or an intense 
beam from muon storage ring l3|. If O13 is not too small, it is likely that leptonic CP 



violation is first explored by LBL experiments with conventional superbeam 13l |. 



To measure CP-violating phase the LBL experiments must run not only with the neutrino 
mode but also with the antineutrino mode. Apart from the problem of parameter degeneracy 
[1^ Esl . Kdi 17, [1^, these measurement would allow us to determine the CP- violating 
phase 5 to a certain accuracy. In the Japan Proton Accelerator Research Complex (J-PARC) 

Hyper-Kamiokande project with upgraded 4MW beam of 50GeV accelerator at J-PARC, 
the accuracy of determination of 6 is expected to be ~ 20 degrees at 3cr CL 0]. 

Running the experiment with antineutrino mode, however, is possible only by overcoming 
a variety of difficulties, much greater ones compared with those in neutrino mode opera- 
tion. Even if we ignore the issue of slightly less intense vr" beam compared to vr"*" beam, 
the antineutrino cross sections are smaller by factor of ~ 3 than neutrino cross sections, 
which results in three-times longer period of data taking, 6 years of u-mode compared to 2 
years of u-mode operation in the J-PARC —>■ Hyper-Kamiokande (hereafter abbreviated as 
JPARC-HK) experiment. Moreover, the background in z/g appearance detection, according 
to the current estimate, are larger by factor of ~ 2 compared with those in z/g detection. 
Hence, antineutrino-mode measurement may be better characterized as an independent ex- 
periment rather than the in-situ measurement. Considering three times longer running time 
it is certainly worthwhile to think about an alternative which can run simultaneously with 
neutrino-mode superbeam appearance experiments. 

In this paper, we point out that a reactor experiment can serve for such purpose. We 
demonstrate that reactor experiments for measuring ^13 with reasonable assumptions on their 
systematic errors can uncover the leptonic CP violation when combined with high-statistics 
neutrino- mode superbeam experiments. In fact, we have pointed out such possibility in our 
previous communication [2^, in which we have demonstrated the complementary role of 
the reactor and the LBL accelerator experiments in determination of the remaining neu- 
trino mixing parameters. The treatment in this paper quantifies our proposal and thereby 
complements and further strengthen our viewpoint of the LBL-reactor complementarity. A 
quantitative treatment of sensitivity for detecting CP violation by reactor-LBL combination 
was also attempted in Ref. jij] but with no indication of signal. See Ref. jl^, 2^1 for 



detailed description of possible experimental designs for reactor experiments for measuring 

^13- 

We remark that the sensitivity to CP violation by our reactor-LBL combined method 
suffers from the problem of parameter degeneracy. However, it can be arranged so that 
it is insensitive to the intrinsic parameter degeneracy 3|. If the superbeam experiment 



is done at the oscillation maximum the combined measurement will allow us to determine 



2 



sin (5. Obviously, the measurement by itself cannot resolve the ambiguity S ^ it — 5, but it 
does not produce a fake CP violation. We have to note that our method suffers from the 
problem of degeneracy due to unknown sign of Am|^ 1^. Even in the case of the JPARC- 
HK experiment in which the matter effect is only modest, it does affects the CP sensitivity 
because the degenerate solutions of 6 differ by ~ 7r/2 in overlapping region of two ellipses 
in the bi-probability plot 3|- Therefore, it is important to know the sign of Amg^ prior to 
the reactor-LBL measurement of S. While the octant ambiguity of ^23 niay also affect 
the CP sensitivity, we do not try to elaborate this point in the present paper. See, however, 
(1) in the concluding remarks. 

We emphasize that reactor experiment cannot replace the antineutrino-mode superbeam 
experiments. It is because the reactor-LBL combined method can detect leptonic CP vi- 
olation only up to ~ 2a CL. Nevertheless, we believe that such reactor-LBL combined 
measurement has a great merit. It will give us the first grip of the structure of leptonic CP 
violation. It will also merit the then ongoing neutrino-mode and the following antineutrino- 
mode superbeam experiments themselves; Even a rough knowledge of the feature of CP 
violation would be very helpful to optimize the setting (such as relative time sharing of u 
and u modes) of the difficult and extremely long-term experiment. 



II. REACTOR-LBL COMBINED MEASUREMENT OF CP VIOLATION 

The principle of detection of leptonic CP violation in a reactor-LBL combined measure- 
ment is very simple. First let us remind the readers the characteristic features of reactor 
measurement of 613. As we have discussed in length in Ref. ^^l, reactor experiment can 
serve for pure measurement of ^13 assuming that Am|]^ is accurately determined by disap- 
pearance measurement of -P(t'^ — > i^fi) in LBL experiments. Namely, it is not contaminated 
by uncertainties due to unknown CP phase 6, the matter effect, and possibly to the octant 
ambiguity 623 11/2 — 623 from which Uf, appearance measurement by LBL experiment 
suffers. 

Now LBL z/g appearance experiment will observe the neutrino oscillation probability 
P{^fi ~^ ^e)- In leading order in /S.m\^/ /S.m\^ it takes the form [2^ 

P{U^ Z/e) = P{y)± = + F±Si3 COS (^5 ± + P0, (1) 

where ± refers to the sign of Amg^^. The coefficients X±, Y±, and Pq are given by 

X.^44,(|^)^.n^(f). (2) 

A2i\ MaA . (aL\ . ( 
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with 



Po = 4,sm'2e,2[^] sin^( — ) (4) 



A,, = ' and P± = A3i±aL, (5) 
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where a = \/2GfNi. denotes the index of refraction in matter with Gp being the Fermi 
constant and Nf, a constant electron number density in the earth. We use in this paper 
the standard notation of the MNS matrix The mass squared difference of neutrinos is 
defined as Am|j = — mf where irii is the mass of the i-ih eigenstate. 

There exist number of reasons for tuning the beam energy to the oscillation maximum 
Ai3 = TT in doing the appearance and the disappearance measurement in LBL experiments, 
as listed in [l^. In this case, cos [5 ± = =f sin5 and ((H) can be solved for sin 5 as 

„„i.£M^ik^i^i£!i. (6) 

We note that, since Oi^ can be measured by reactor experiments, the right-hand side (RHS) 
of © consists solely of experimentally measurable quantities. Therefore, LBL measurement 
of Pi^fj, —>■ Ve)-, when combined with the reactor experiment, implies measurement of sin 5. 

In the rest of this paper, we try to elaborate our treatment by including suitably estimated 
experimental uncertainties of both LBL and the reactor experiments. As indicated in 0, 
the accuracy of measurement of sin 5 solely depends upon how precisely -P(z^) and sis in 
the RHS can be determined in LBL and reactor experiments, respectively. We take the 
best possible case among the concrete proposals of LBL experiments currently available in 
the community, the JPARC-HK experiment assuming 4MW beam power and 540 kton as 
the fiducial volume of the detector jl^]. However, most probably our conclusion does not 
heavily depend on any detailed experimental setting in the particular experiment, once the 
accuracy of measurement of the Uf. appearance probability reaches to that level and if the 
baseline is not too long. For the reactor experiment, we present our results in units of 
GWth ■ ton • year exposure to allow application to wider class of experiments. Our results 
may be useful to indicate what condition must be met to uncover the leptonic CP violation 
in such reactor-LBL combined measurement. 



III. TREATMENT OF ERRORS IN LBL AND REACTOR EXPERIMENTS 

To carry out quantitative analyses of the sensitivity for detecting CP violation, we must 
first establish the method for statistical treatment of LBL and reactor experiments. 

A. Treatment of errors in the JPARC-HK experiment 

We consider neutrino-mode appearance measurement for 2 years in the JPARC-HK ex- 
periment. For definiteness, we use the neutrino fiux estimated for the off-axis 2° beam j^. 
We define Ax^ for the experiment as 

Xj.PARC. - ^best + jVbg + cri(Arb-t)2 + al^{N^G? ' ^ ^ 

where and A'bg represent the expected number of signal and background events, respec- 
tively, computed with the cross section in [i^. A^^*''^* is defined as the number of signal 
event Ni, calculated with the best-fit values of the "experimental data", which is to be 
tested against the CP conserving hypothesis, 5 = 0. as and (Tbg represent the fractional 
uncertainties of the estimation of the number of signal and background events, respectively. 
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between detectors 


single detector 


correlated 


uncorrelated 


between bins 


correlated 


(JDB = 2.5% 


cJdB = 0.5% 


CTB ~ 2.6% 


uncorrelated 


cTDb = 2.5% 


(jdb = 0.5% 


cJb ~ 2.6% 


total number of events 


cJD ~ 2.6% 


CTd ~ 0.5% 


dsys ^ 2.7% 



TABLE I: Listed are assumed values of systematic errors cjdBi CDb) crdB, and cJdb. The subscripts 
D (d) and B (b) are represent the correlated (uncorrelated) error among detectors and bins, re- 
spectively. Using those four values, the errors for the total number of events and for single detector 
are calculated. 

Following j^, we use as = ctbg = 2% in our analysis. (See Sect. HYlfor more about how to 
use Ax^ in our procedure to determine the sensitivity region for CP violation.) 

While we do not use the spectrum information in a direct way in our analysis, we need to 
estimate how the experimental event selection affects the spectrum to calculate the number 
of signal and background events. The most important cut is to suppress the background 
events due to 7r°. We use the simulated spectra after the cut calculated by the JPARC-SK 
group [io'l and evaluate the reduction rate due to cut in each energy bin of 50MeV width. 
The procedure is applied to calculate the number of signal after the cut for any values of 
mixing parameters. In this way, the total numbers of events within energy range 0.4-1.2GeV 
are calculated and used in our analysis. 



B. Treatment of errors in the reactor experiment 

In this paper we consider the case of single reactor and two (near and far) detector 
complex.^ The far detector is placed 1.7km away from the reactor, the optimal distance for 
lAm^^l = 2.5 X lO^^eVl We assume that a near detector identical with the far detector is 
placed 300m away from the reactor to reduce systematic errors.^ 

We consider four types of systematic errors: ctdb, crob, c"dB, and (Tdb. The subscript D 
(d) represents the fact that the error is correlated (uncorrelated) between detectors. The 
subscript B (b) represents that the error is correlated (uncorrelated) among bins. To indicate 
nature of these respective errors, we list below some examples of the errors in each category: 

(Tdb: error in estimation of reactor power 

a-Di,: error in estimation of detection cross sections 

a^B'- error in estimation of fiducial volume of each detector 

(Tdb: errors inherent to detectors such as artificial firing of photomultiplier tubes 
Although the values of a^B for far and near detectors, for example, can be different from 
each other, we neglect such difference for simplicity. The values of systematic errors we 
assume are listed in Table HI 



^ The current proposal by the Japanese group 12.^ 1 plans to utilize flux from 7 reactors observe(i by two 
near and a far detectors. It is shown even in this case that an effective 1 reactor-2 detector approximation 
gives a very good estimation of the sensitivity |3ol |. 

^ The closer the near detector to reactor, the better the sensitivity in the single-reactor case because of 
larger number of events. The situation is, however, more subtle for multiple-reactor case [3o| . 
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As will be briefly explained in the appendix, the errors ctd and for the total number 
of events are obtained as 

where a = n, f are the index for near and far detectors, and i runs over number of bins. We 
use 14 bins of O.SMeV width in l-8MeV window of visible energy, E^isi = Ep^ — O.SMeV. 
The coefficient of and a^^ is about 1/9 in our analysis almost independently of a. Since 
relative normalization errors are \/2 times of uncorrelated errors, cxd — 0.5% is consistent 
with the value used in ji^. In Ref. [i^, the most pessimistic assumption ctdb = c^dB = 
was taken for bin-by-bin distribution of errors. The value of (jgyg = a^ + a^ is also consistent 
with the total systematic error of the CHOOZ experiment. In summary, we feel that the 
errors listed in Table HI are not too optimistic ones and are likely to be realized in the setting 
discussed in HI 113,13. 

Our definition of Ax^^^^^^. is 



^xLct = min 



react 

^ [ {Ng, - (1 + + + a)iVbest)^ ^2 

2.1 ^w + ^2^(^best)2 + a 



as 

a=f,n 



Db I "^dB 



where Nai represents the theoretical number of events at a-detector within i-th bin. Again, 
N^l'^^ is defined as the number of signal event calculated with the best-fit parameters of the 
"experimental data" . The minimization in is achieved analytically, and then we obtain 

AxLct = (^^^y>-^(^l), (10) 



X 



>/i-A^r ^ ^T_( Nm-NT , .... 

' y — Tvrbest ' ' \^^/ 



Arbest ' ' AThe 



est ' 



+ ^i/28 + al^ ^ J + al^H^s, (12) 

where /„ represents the n x n identity matrix and Hn represents the n x n matrix whose 
elements are all unity. Notice that an infinitely good sensitivity is obtained for infinite 
number of events if cxdb vanishes because det(VQ_goes to zero for the case which explains the 
apparently curious behavior seen in Fig. 2 of '21]. See Q for more about the equivalence 
between the "pull" and the covariance matrix methods. 

To indicate the expected sensitivity of the reactor experiment with the systematic errors 
listed in Table HI we present in Fig. [T]the excluded region in sin^ 2^i3-| Amg^^l space in the 
absence of fiux depletion (6'5'|^* = 0) for 10^, 4 x 10^, and 10^ GWth ■ ton ■ year exposure. 
The z/g detection efficiency of 70% is assumed j3,|2^. The number of events expected during 
these exposure are about 10^, 4 x 10^, 10^ Ue events, respectively, at the far detector.^ Notice 



V = diae 



^ In the rate-only analysis without binning, the sensitivity is saturated at the number of i/^ events around 
10^ 
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that what we mean by numbers in units of GWth is the thermal power actually generated 
from reactors and it should not be confused with the maximal thermal power of reactors. 
Assuming average 80% operation efficiency the above three cases correspond approximately 
to 0.5, 2, and 5 years running, respectively, for 100 ton detector at the Kashiwazaki-Kariwa 
nuclear power plant whose maximal thermal power is 24.3GWth- 

IV. ESTIMATION OF SENSITIVITY OF REACTOR-LBL COMBINED DETEC- 
TION OF CP VIOLATION 

To estimate the sensitivity of the reactor-LBL combined measurement to leptonic CP 
violation, we define the combined Ax^ as 

Axlpi{S;S^^'\sm^2e^f) = min Axcp{S,sm^ 2613] S^'''\sm^ 26^^) 

sin^ 2^13 

^ min {AxlpARc.(<5,sin2 2^i3;5''=^\sin2 2^^f) 

+ AxLct(sm^2^i3;sin22^^r)}. (13) 

We take the following procedure in our analysis. We pick up a point in the two-dimensional 
parameter space spanned by 5^*^^* and sin^ 29\^^^ and make the hypothesis test on whether 
the point is consistent with CP conservation within 90% CL. For this purpose, we use 
the projected Ax^ onto one-dimensional 6 space, Axcpi, as defined in (|13|) and then the 
statistical criterion for 90% CL is Ax^p^^ < 2.7. Then, a collection of points in the parameter 
space which are consistent with CP conservation form a region surrounded by a contour in 
^best_gjj^2 2^best gpg^^g^ g^g -^jn Iqq showu in Figs. 121121 below. 

The neutrino mixing parameters are taken as follows: |Am|^| = 2.5 x 10~'^eV^, Am2i = 
7.3 X 10-^eV^ tan2 6'i2 = 0.38, and sin^ 26*23 = 1- Notice that the high-Am^i LMA-II solar 
neutrino solution is now excluded at 3a CL by the global analysis of all data with reanalyzed 
day-night variation of fiux at Super-Kamiokande [32|, and at 99% CL by the one with SNO 
salt phase data [3^. The earth matter density is taken to be p = 2.3 g ■ cm~^ and the 
electron number density is computed with electron fraction = 0.5. 

A. CP sensitivity in the case of known sign of Am^^ 

In Fig. 121 the regions consistent with CP conservation at 90% CL are drawn for Arrig^ > 
case in the region — 7r/2 < 6^^^^ < 7r/2. The thin-solid, the solid, and the thick-solid lines 
are for 10^, 4 x 10^, and lO^GWth ■ ton ■ years, respectively, and the regions consistent with 
CP conservation are within the envelope of these contours.'* We remark that the present 



Since we rely on hypothesis test with 1 degree of freedom (1 d.o.f.) the information of sin^ 26'i3 is lost 
through the process of minimization in H13I) . The individual contours presented in Fig. [21 indicate the 
region Axcp < 2.7 for eight assumed values of sin^ 26*13 which range from 0.02 to 0.16. In this way the 
figure is designed so that the envelop of the contours gives the region of CP conservation at 90% CL 
by 1 d.o.f. analysis, and at the same time carries some informations of how the sensitivity regions are 
determined by the interplay between the reactor and the LBL measurement. We hope that no confusion 
arise. 
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constraint on ^13 becomes milder to sin^2^i3 < 0.25 at Scr CL js^ by the smaller values 
of |Am|j^| indicated by the reanalysis of atmospheric neutrino data [365. Notice that the 
other half region of S^^^^ gives the identical contours apart from tiny difference which arises 
because the peak energy of the off-axis 2° beam is slightly off the oscillation maximum. 

If an experimental best fit point falls into outside the envelope of those regions, it gives an 
indication for leptonic CP violation because it is inconsistent with the hypothesis S = at 
90% CL. We observe from Fig.|21that there is a chance for reactor-LBL combined experiment 
of seeing an indication of CP violation for relatively large e^l''\ sin^ 26^1''^ > 0.03 at 90% CL. 
We believe that this is the first time that a possibility is raised for detecting leptonic CP 
violation based on a quantitative treatment of experimental errors by a method different 
from the conventional one of comparing neutrino and antineutrino appearance measurement 
in LBL experiments. 

The sign of Amg^^ is taken to be positive in Fig. |21 which corresponds to the normal mass 
hierarchy. If we flip the sign of Amgj^ (the case of inverted mass hierarchy) we obtain almost 
identical CP sensitivity contours. It is demonstrated in Figs. 12^^ and|2b which serve also for 
the discussion in the next subsection. By comparing the contours depicted by thick-solid 
and thick-dashed lines in Figs. (Am|^ > 0) andEb (Amg^^ < 0), respectively, it is clear 
that the CP sensitivity is almost identical between positive and negative Aml^. The largest 
noticeable changes are shifts of the end points of the contours toward smaller (larger) S in 
the first (fourth) quadrants by about 10% (a few %) at sin^ 2^13 = 0.1. Namely, the both 
end points slightly move toward better sensitivities for the inverted mass hierarchy. 

The sensitivity contour of CP violation is determined as an interplay between constraints 
from reactor and accelerator experiments. The former gives a rectangular box in the 6^'^^^- 
sin^ 29\^^^ space, whereas the latter gives the equal-P(i/) contour determined by ^ under 
the hypothesis 6 = with finite width due to errors, as indicated in Figs. El and El In region 
of parameter space where both of these two constraints are satisfied, the best fit parameter 
is consistent with CP conservation. Outside the region the CP symmetry is violated at 
90% CL. The discovery potential for CP violation diminishes at small sin^ 29i^^^ primarily 
because -P(z^) becomes less sensitive to 6 at smaller 6'i3, while the reactor constraint on 
sin^ 26'5^|*** is roughly independent of 613 |20 |. 

B. CP sensitivity in the case of unknown sign of Amg^^ 

So far we have assumed that we know the sign of Am|^ prior to the search for CP 
violation by the reactor and the JPARC-HK experiments. But, it may not be the case 
unless LBL experiments with sufficiently long baseline start to operate in a timely fashion. 
In this subsection we assume the pessimistic situation of unknown sign of Amg^^ and try to 
clarify the influence of our ignorance of the sign on the detectability of CP violation by our 
method. 

If the sign of Amg^^ is not known, the procedure of obtaining the sensitivity region for 
detecting CP violation has to be altered. It is because we have to allow such possibility as 
that we fit the data by using wrong assumption for the sign. In Fig. (jSt") we present 
the results of the similar sensitivity analysis for detecting CP violation as we did in the 
previous section by assuming that the sign of (Amg^^)^®*^*, which is chosen by nature, is 
positive (negative). It is obvious from Fig. ^ (jSb) that the contours of CP conservation 
moves to rightward (leftward) if the wrong sign is assumed in the hypothesis test, essentially 
wiping out about half of the CP sensitive region of 5^*^'^*. 
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The results can be confusing and some of the readers might have naively interpreted, 
by combining Figs. and Eb, that there is no sensitivity region in S^"^^^- sin^ 2^2^* plane. 
To resolve the puzzling feature we present in Fig. |3] the regions which are consistent with 
CP conservation by contours in the plane spanned by observable quantities, the number of 
events in the reactor and the JPARC-HK z/g appearance experiments.^ This plot indicates 
that the sensitivity region for detecting CP violation does not disappear but becomes about 
half. Which region of S is CP sensitive depend upon the sign of Am^g, or in other word on 
the location in bi- number of event plane in Fig. |3J For complete clarity, we have placed three 
different symbols in Fig. |3] and at the same time in Fig. El to indicate which points in the 
space of observables correspond to which points in the CP sensitivity plot. Note that the 
point indicated by a cross in Fig. 0] corresponds to two values of 6^^^^ because of unknown 
sign of Am^g. 



V. CONCLUDING REMARKS 

In this paper we have pointed out a new method for detecting leptonic CP violation by 
combining reactor measurement of ^13 with high-statistics z/g appearance measurement in 
LBL accelerator experiments. A salient feature of our method are that one can perform the 
measurement prior to the lengthy antineutrino running in LBL experiments. We conclude 
with several remarks: 

(1) If 623 is not maximal the parameter degeneracy due to the octant ambiguity of 623 will 
also affect the sensitivity of detecting CP violation. On the other hand, we have discussed in 
our previous communication that the octant degeneracy may be resolved by combining 
reactor measurement of ^13 with the LBL appearance measurement in both neutrino and 
antineutrino channels. It would be very interesting to reexamine the possibility in the context 
of this work to clarify to what extent it cures the further uncertainty in the sensitivity of 
detecting CP violation mentioned above. 

(2) We have examined a pessimistic scenario to run the JPARC-SK experiment with 0.75MW 
proton beam power and the fiducial volume of 22.5kt, while waiting for the construction of 
Hyper-Kamiokande. As is shown in Fig. the sensitivity to CP violation becomes worse 
but still remains for its 10 years running. 

(3) The reactor experiment described in this paper may be regarded as the phase II of the 



currently proposed reactor experiments for measuring ^13 [2^ |23|, |2J| , and how to improve the 
systematic errors should be carefully investigated during running the phase I experiments. 
If it is possible to significantly improve the systematic errors over those given in Table Ql it 
may be possible to extend the CP sensitivity to the region sin^ 2^13 < 0.03. 

(4) From Figs. 121121 it is likely that detection of CP violation requires ~ 10'^ GWth ■ ton ■ year 
measurement by the reactor experiment. Now there is a choice between two options: stronger 
power source with smaller detectors, or weaker power source with larger detectors. If these 
is no natural or existing holes with enough overburden for the detectors the first option 



^ Notice, however, that we have used binned data, not merely the total number of events, in analyzing 
reactor experiment to obtain the contours. 
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might be more advantageous because larger detectors require deeper hole to keep the signal 
to noise ratio equal. 
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APPENDIX A: 
COMPARISON 



CANCELLATION OF ERRORS BY NEAR-FAR DETECTOR 



This appendix is meant to be a pedagogical note in which we try to clarify the feature 
of cancellation of systematic errors by near-far detector comparison and the relationship 
between over-all and bin-by-bin errors. 

The definition of Ax^ for two detector system is 



nf 



min Ax,^ 



nf 



a] 



mm 



[Nf - (1 + a)A^}^^*}^ 



|iV„-(l + a)iVbe^t|2 



+ 



a 
+ — 



(Al) 



where Nf [Nn) is the theoretical total number of events expected to be measured at far (near) 
detector. The quantities with superscript "best" are defined as the ones calculated with the 
best-fit values of the "experimental data" , which are to be tested against the CP conserving 
case. (Td and cxd are correlated and uncorrelated errors between detectors, respectively. 

We discuss statistical average of an observable O by the Gaussian probability distribution 
function as 



<0>i 



C j dNfdNndaOexp (^-^Axlj 



.a 



(A2) 



where C is the normalization constant to make < 1 > unity. Note that the integration with 
respect to a is equivalent to the minimization in ()A1|) . After the minimization, it takes the 
following form which is generic to the Gaussian distribution. 



< X > < xy > 

< yx > < y"^ > 



(A3) 



X 



Nf - A^best 
jybest • 



Nn-N. 



best 



jybest 



(A4) 
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In order to examine the feature of near-far cancellation of errors, it is valuable to transform 



X and y as 



/ n f n 



Then, ^Xnj can be written as in the form ()A3|) with 

<^'> = i + i + 2c^d, (A6) 

/ " 

<Y'> = ]^ + ]^ + 2a^ + 44, (A7) 
/ " 

<XY> = <YX>= — ^ - — (A8) 

/ n 

It is evident in ()A6|) that the correlated systematic errors cancel by the near-far comparison. 
The systematic error V2ad in < > is referred in |2ol l37l | as the relative normalization 
error. 

We briefly treat the case of two bins with infinite statistics to illustrate the importance 

(.nf 

-1 



of uncorrelated errors. In this case X subspace of ^Xnf can be written as 



It is clear that a^h = leads to the diverge of ^Xnj except for the best fit point (Xj = Yi = 0), 
which means that the infinite precision can be achieved for the case. Thus, adb must be 
treated with great care. 

Next we derive the relationship between over-all and bin-by-bin errors that was used in 
the text, (jHl). For simplicity, we consider the case of one detector with two bins. Then, Ax^ 
for the case is defined as 

Axi2 = minAx?2(«) 

(AlO) 



mm 

a 



{Ni - {1 + a)N]''''} {^2 - (l + a)X2''*} 



where Ni and N2 are the expected numbers of events within first and second bins, respec- 
tively, and (Tb {<^h) denotes the correlated (uncorrelated) error between bins. 
To obtain the error for the total number of events, we define 

, AT _ Afhest ^ 

-tot-E^^^' AT--. 5^ AT-- (All) 



'tot 



Then, we obtain 



2 , 

^best- B+ b ^^bcst^2 



+ + K ^Z^LJ • (A12) 

tot 
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One can show that the same treatment goes though for arbitrary number of bins. The 
coefficient of cr^ is almost 1/9 in our analysis (14 bins). 
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FIG. 1: The excluded regions at 90% CL in the absence of v^, disappearance (^xf^* = 0) are 
drawn for 10^, 4 x 10^, and 10^ GWth • ton • year exposure of a reactor experiment by thin-sohd 
(red), solid (green), and thick-solid (blue) lines, respectively. The far (near) detector is placed 
1.7km (300m) away from the reactor. We assume that |Am3;^| is precisely measured by LBL 
experiments and adopt the analysis with one degree of freedom (Axreact = ^ •''')• ^^e the value 
lAmgil = 2.5 X 10^'^eV^ as indicated by the dashed-doted line in the figure. In our analysis, we 
use 14 bins of 0.5MeV width in l-8MeV window of visible energy with the systematic errors listed 
in Table ni 
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FIG. 2: The contours are plotted for eight assumed values of sin^ 2^13 which range from 0.02 to 
0.16 to indicate the regions consistent with the hypothesis 5 = at 90% CL (Ax^p = 2.7) by 
the reactor-LBL combined measurement. If an experimental best fit point falls into outside the 
envelope of those regions, it gives an evidence for leptonic CP violation at 90% CL. The thin-solid 
(red), solid (green), and thick-solid (blue) lines are for 10"^, 4x 10"^, and lO^GWth -ton-year exposure 
of a reactor experiment, respectively, corresponding to about 0.5, 2, and 5 years exposure of 100 ton 
detectors at the Kashiwazaki-Kariwa nuclear power plant. For the JPARC-HK experiment, 2 years 
measurement with off-axis 2° z/^ beam is assumed. (See the text for more details.) The normal 
mass hierarchy, Am|-^ > 0, is assumed. 
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FIG. 3: The contours which surround the region consistent with CP conservation are plotted in 
Fig. (jnb) by assuming (Am^i)^^^^ > {{Am'^i)^'^'^^ < 0) as nature's choice. If the right (wrong) 
sign is used as the hypothesis with 5 = 0, the contours indicated by the thick (thin) hues result in 
both figures. The three symbols, a cross, open and solid circles are placed on the figures as well 
as in Fig. to indicate the relationship between observed numbers of events and the results of CP 
sensitivity analysis. 
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Nbest/4MW-540kt-2yr 

FIG. 4: The contours which surround the region consistent with CP conservation are plotted in 
the number-of-events space of reactor and the LBL experiments. The thin hues correspond to 
jbest _ ji. The three symbols, a cross, open and solid circles are placed on the figures as well 
as in Fig. to indicate the relationship between observed numbers of events and the results of CP 
sensitivity analysis. 
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FIG. 5: The same as in Fig. ^ but with measurement in 10 years running of the JPARC-SK 
experiment with 0.75MW beam power and 22.5kt detector (Super-Kamiokande). Although each 
contour becomes thicker because of a factor of ~ 25 lower statistics of the experiment, the sensitivity 
to CP violation still exists at 90% CL. 
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